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Possible signature of hypernova nucleosynthesis in a 
beryllium rich halo dwarP 
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ABSTRACT 

As part of a large survey of halo and thick disc stars, we found one halo star, HD 
106038, exceptionally overabundant in beryllium. In spite of its low metallicity, [Fe/H] 
= —1.26, the star has log(Be/H) = —10.60, which is similar to the solar meteoritic 
abundance, log(Be/H) — —10.58. This abundance is more than ten times higher the 
abundance of stars with similar metallicity and cannot be explained by models of chem- 
ical evolution of the Galaxy that include the standard theory of cosmic-ray spallation. 
No other halo star exhibiting such a beryllium overabundance is known. In addition, 
overabundances of Li, Si, Ni, Y, and Ba are also observed. We suggest that all these 
chemical peculiarities, but the Ba abundance, can be simultaneously explained if the 
star was formed in the vicinity of a hypernova. 



Key words: 

106038 



stars: abundances - stars: chemically peculiar - stars: individual: HD 



1 INTRODUCTION 

' The single stable isotope of beryllium, ^Be, is a pure prod- 
uct of cosmic-ray spallation of heavy (mostly CNO) nuclei 
eeves. Fowler fc HoyleHigTOj). Analyses of Be abundances 
in m etal poor stars ( Molaro et al.l Il997l : iBoesgaard et al.) 
1 19991 ) have found a relationship between [Fe/Hj3 and 
log(Be/H) with slope close to one, and between [0/H] and 
log(Be/H) with slope between 1 and 1.5, depending on the 
oxygen indicator used. Independently of the behaviour of 
the [O/Fe] ratio at lower metallicities, these re sults sugges t 
a primary production of Be in the early Galaxy l|Kinell200ll ). 

As a primary element, and assuming cosmic-rays to be 
globally transported across the early Galaxy, Be may show 
a s maller scatter than the products of stellar nucleosynthe- 
sis l|Suzuki fc Yoshiill200ll ) at a given time, suggesting its 
potential use as a cosmochronometer. 



So far, the linear relations appear to be very tight, show- 
ing a surprisingly low scatter comparable to the measure- 
ment errors. This picture, however, might change with the 
increase of the samples anal ysed, as hinted by the results of 
IBoesgaard fc Novickil (|2006l ). Nevertheless, turn-off stars of 
the globular clusters NGC 6397 and NGC 6752 were found 
l|Pasquini et al.ll2004 [20071 ) to have the same Be abundance 
of field stars of the same metallicity. This strongly sup- 
port the production of Be to be a global process. Ages de- 
rived from these abundances, in a comparison with a model 
of the evolution of Be with time l|Valle et al.ll2003 ). show 
an excellent agreement with ages derived from theoretical 
isochrones, su pporting the us e of B e as a cosmochronome- 
ter. Moreover, IPasquini et al.l l|2005h showed that Be abun- 
dances could be used to study the differences in the time 
scales of star formation in the halo and the thick disc of the 
Galaxy. 



* Based on observations made with ESO VLT, at Paranal Ob- 
servatory, under programmes 076.B-0133 and on data obtained 
from the ESO/ST-ECF Science Archive Facility 
t E-mail: rodolfo@astro.iag.usp.br (RS); lpasquin@eso.org (LP) 
1 [A/B] = log [N(A)/N(B)], - log [N(A)/N(B)]0 



In this letter, we report the discovery of an extremely 
Be enriched halo star, HD 106038, with an abundance 1.2 
dex higher than stars of similar metallicity. This unique star 
deviates considerably from the observed relations of Be with 
Fe and O. It was identified during the analysis of a large 
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Figure 1. Comparison between the spectra of HD 106038 (solid 
line) and of HIP 7459 (dashed line), a star with elose atmospheric 
parameters and similar metallicity, in the Be region. The domi- 
nating element of the nearby blended features are also indicated. 
The V and Ti features have the same strength in the two stars 
while some difference in the Zr line is noted. 



Figure 2. The abundance of Be as a function of [Fe/H], The 
stars HD 106038 and HIP 7459 (filled squares ) are compared 
to the linear relation defined by the stars from iBoesgaard et al.l 
ll999l) (open circles). The starred symbols are t he st ars from 
iBoeseaard et ah 1 lll999l ') and lBoeseaard fc Novick i 1 I2OO6I) . Two of 
them deviate from the linear relation by ~ 0.50 dex. 



sample containing near to one hundred halo and thick disc 
stars (Smiljanic et al. 2008, in preparation). 

Neither the standard scenario for Be production, in- 
volving spallatio n of cosmic-rays on nuclei of the inter- 
stellar medium llValle et al.l |2002| ). nor the superbubbles 
(SBs) scenario 



ParizotI bOOoF seem to be able to pro- 



duce such Be enriched object s. The SBs model pre dict a 
scatter in the Be abundance ("Pari zot fc Drurvll200(il) that 
may explain the stars found by iBoesgaard et al. ( 1999 ) and 
IBoesgaard &: Novickil (|2006l ) which have similar atmospheric 
parameters but Be abundances differing by ~ 0.5 dex. The 
very high Be abundance in HD 106038, however, would re- 
quire an extremely poor mixing of the SN e ejecta with the 
ISM which seems to be difficuH to justify (|Parizotll2000l ). 

To the best of our knowledge, there is only one other 
case of extremely Be enhanced star in the literature The 
star J37 of the op en cluster NGC 6633 was found by 
lAshweU et al] l|2005l ) to have log (Be/H) = -9.0 ± 0.5. The 
chemical peculiarities of star J37 might be best explained 
by the a ccretion of rocky mat erial similar to chondritic me- 
teorites (|Ashwell et al.l 20051 ) . As we shall see below, the 
accretion of such a material is unlikely for our population II 
star. 



2 DATA AND ANALYSIS 

The science raw data and calibration files of HD 106038 were 
retrieved from the ESO science archive facility. The spectra 
were originally obtained in 2000 April 12 with the UVES 



^ We exclude from the discussion the chemically peculiar A or F 
stars with enhanced Be lines. The peculiar abundances of these 
stars are thought to be caused by e ffects of difi^usion. As shown by 
iRichard. Michaud fc Richer] 1 I2OO2I) . these eff'ects do not result in 
overabundances in stars with similar temperature and metallicity 
as HD 106038. 



spectrograph (jPekker et al.ll2000l ) of the ESO VLT at Cerro 
Paranal, Chile. The data of HIP 7459 (CD-61 282), a halo 
star used as comparison in this work, were obtained in 2005 
September 22 with the same instrument. The spectra have 
R ~ 45000 and a final S/N ~ 70 in the Be region. 

For b oth stars, we adopt the a tmospheric parameters 
derived by iNissen fc Schusted l|l997h . The parameters were 
calculated with the standard spectroscopic approach using 
Fe I and Fe II lines (Table [TJ . We refer the reader to the 
original work for more details. 

Abundances (Table [l| were derived through the syn- 
thesis of the spectrum around the Be II resonance lines 
at 3131.065 A and 3130.420 A. The codes for calculat - 
ing synthetic spectra are described in ICoelho et all l|2005h . 
Grids of model atmospheres wit hout overshooting calcu- 
lated by the ATLAS9 program (|Castelh fc Kurud I2OO3I ) 
were adopted . The list of atomic lines is that compiled by 
iPrimas et al.l Il 997) and the molecular line list is described 



Coelho et al.l (|2005D . A solar Be abundance was derived 



using the UVES solar spectrum. We estimate the total un- 
certainty from atmospheric parameters, continuum place- 
ment, and synthetic fit affecting the Be abundance to be a 
— ± 0.13 dex. Zero point errors might also be present due, 
for example, to the adopted model atmosphere. However, 
we are conducting a differential analysis and these should 
cancel out in a comparison between similar stars. 

A comparison between the spectra of the two stars is 
shown in Figure [T] in which the extreme enhancement of the 
Be lines of HD 106038 when compared to the normal HIP 
7459 is clear, confirming its very high Be overabundance. 

We show in Figures[2]and|3]these two stars in the [Fe/H] 
vs. log (Be/H) and [0/H] vs. log ( Be/H) diagrams, r espec - 
tive lly, together with the stars fr om IBoesgaard et al.l (Il999l ) 
and IBoesgaard fc Novicki (I2OO6I). Of pa r ticula r interest are 
the two stars from IBoesgaard fc Novickil (|2006l ) that deviate 
from the linear trend. The anomalous position of HD 106038 
clearly stands out. 
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Figure 3. The abundance o f Be as a function of [O/H]. Symbols 
are as in Fig. [2] The stars of lBoeseaard fc Novickil (120061) are not 
shown since O abundances were not derived by these authors. 



Table 1. The adopted atmospheric parameters and beryllium 
abundances derived using synthetic spectra for HD 106038, HIP 
7459, a comparison star, and for the Sun. 



star 


Tcff 


log g 




[Fc/H] 


log(Be/H) 




K 




km s^^ 






Sun 


5777 


4.44 


1.00 


0.00 


-10.9 


HD 106038 


6046 


4.46 


1.34 


-1.26 


-10.6 


HIP 7459 


5909 


4.46 


1.23 


-1.15 


-11.9 



2.1 Chemical abundances in the Hterature 

Information on other chemical abundances might help un- 
derstanding the origin of the Be overabundance. An over- 
abundance of CNO elements, for example, would offer a good 
explanation for the enhancement, since these elements are 
do minant in the producti on of Be by spallation processes. 

lAsplund et al.l (|2006l ) determined a lithium abundance 
of A(''Li) = 2.48 and claimed a detection of ''Li, ''Li/^Li = 
0.031, compatible with the other 8 detections out of a sample 
of 26 stars. The high ^Li abundance, however, results in a 
high ''Li abundance, A(''Li) = 0.97, while the mean for the 
other detections is 0.80. 

Its ^Li is par ticularly re markable since it is larger than 
th e Spite plateau (ISpite fc Sp ite 1982). The plateau as found 
bv lAsplund et al.l l|2006l ') is at 2.22, which means HD 106038 
has an abundance excess of A^(^Li) = 2.13. Given that 
lithium is also expected to be produced by cosmic-ray spal- 
lation it seems safe to conclude that both '^Li and Be over- 
abundances have the same origin. 

Abundances of other elements available in the litera- 
ture are shown in Figure |4l Abu ndances of Na, Ca, Mg, S i, 
Ti, Cr, Ni, Y, and Ba are from iNissen fc Schuste3 lll997^ . 
The carbon abu ndance is taken from Fabbian et al.l ( 20061 ) 
and oxygen from iMelendez et" "ail ll2q06l ). Abundances of Mn 
and Sc are fr om Nissen et al. (|2000l ) and of S and Zn from 
iNissen et al.l (|2007l ). In the same figure we also show the 
mean abundances of the samples analysed in these same 
works for comparison. 



Figure 4. Elemental abundances, [X/Fe], of HD 106038 deter- 
mined in the literature. The abundances are represented by the 
full squares. In this case the error bar denotes the actual uncer- 
tainty quoted by the original work. The starred symbols indicate 
the mean abundance of that element for the remaining sample as 
found in the same work. In this case the error bar indicates the 
range of abundances for that given element in the original work. 
HD 106038 is clearly overabundant in Si, Ni, Y, and Ba, and show 
slightly larger abundances of C, S, Mg, and Zn when compared 
to the mean of the original samples. 



In addition to Be and Li, the star also shows clear en- 
hanced abundances of Si, Ni, and of the neutron capture el- 
ements Y and Ba. An enrichment in s-process elements may 
explain the enhanced Zr line in Figure [1] We also note the 
possibly larger amounts of C, S, Mg, and Zn, though these 
remain marginally compatible with the mean abundances of 
the samples. 



3 DISCUSSION 

Since the standard scenario for cosmic-ray spallation does 
not explain the enhancement of Be in HD 106038, a pecu- 
liar and/or rare event may be related to its formation. A 
combination of two or more rare events to produce the ob- 
served features is unlikely, we therefore concentrate on single 
events. 

To reproduce the very particular chemical pattern of 
HD 106038, a nucleosynthetic site must be able to over- 
produce Si and Ni without overproducing other a and iron 
group elements. Elements in normal halo stars with the same 
metallicity as HD106038 come mostly from SNe II. It is 
therefore unlikely that the same SNe II may produce the 
observed enhanced [Si/Fe] and [Ni/Fe] ratios. Moreover, it 
has about 16 times more Be than what models involving SN 
II predict for its metallicity (jValle et al.ll2002l ). SNe la pro- 
duce large amounts of Fe, Ni, and other iron group elements 
but are not expected to produce large amounts of O and 
Si ijlwamoto et al.l[l999h . Therefore, in case of a contribu- 
tion from SNe la, ratios such as [0/Fe] and [Si/Fe] would 
fall below normal halo stars, contrary to what is observed. 
Moreover, SNe la are expected to produce about one order of 
magn itude less spallation products than SN II |Fields et al.l 
l2002h . 
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All these features, however, may possibly be found in 
the ejecta of a hypernova (HNe), which can be enriched in 
both intermedia te mass ele ments (as S and Si) a s in iron 
group elements (jNakamura et al ., 2001.; Podsiadlowski et"al] 
|2002| ). Moreover, HNe can produce large amounts of Be and 
Li by spallation ([Fields et al.ll200^ : INakamura fc Shigevamal 
12004 ). Therefore, we suggest that the material which formed 
this star was probably contaminated by the nucleosynthetic 
products of a HNe. 



3.1 Hypernova 

Hypernovae are core-collapse SNe (usually of type Ic) with 
exceptionally large kinetic energy production, resulting in 
spectra dominated by very broad absorption line blends 
(|MazzaU et al.ll2000l) . The energy released in the explosion 
can be one order of magnitude larger t han that of normal 
core-collapse SNe ([iwamoto et al.ll2003l ). Some hypernovae, 
typically the most mas sive and energetic ev ents, are linked 

to Gamma-Ray bursts (llwamoto et al.l 19981). 

'Fie lds et al.l (12002^ and'Nakamu ra fc Shigevamal ([20041 ) 
calculated the yields o f spallation products res u lting from 
HNe explosions. While INakamura fc Shigevamal l|2004l ) cal- 
culate the energy distribution of the ejecta with a hydro- 
dynamic code and solve the cosmic-ray transfer equation, 
[Fields et all ([20021 ) use an empirical formula for the energy 
distribution and do not solve the transfer equation but adopt 
an approximation to have the mass fraction of the ejecta that 
produ ces the spallation products. 'Nakamura fc Shigev amal 
20041) claim the simplifications adopted by Fields et al.l 
20021 ) to overestimate the yields by a factor ~ 3. 



The yield of Be per HNe can be one o r two orders 
of m a gnitude larger than the one pe r SNe H ([Fields et al.l 
|2002| ; INakamura fc Shigevamal |2004| ). However, as a rare 
event, they a re not major contributors of Be in the Galaxy. 
iFields et all (|2002l ) predict ^Li/^Be 8.6 and Be/O ~ 
5.6 X 10~^ (both ratios by nu mber). The calculations by 
INakamura fc Shigevamal ([20041 ) predict '^Li/^Be 4.2, also 
by number. Both predictions are close to what is observed in 
HD 106038; the ratio between the observed excess of '^Li and 
^Be is '^Li/^Be = 5.6 (while the HNe is expected to have pro- 
duced all the observed Be abundance, it is responsible only 
for the excess of Li with respect to the primordial plateau). 
We cannot estimate the contribution of the possible HNe on 
the observed oxygen abundance, thus only a lower limit can 
be placed, Be/O >1.9 x lO"'^, given by the assumption that 
all the observed oxygen has been produced by the HNe. 

Both models, howev er, predict much m ore ®Li than oh- 
served ^Li/^Li 1.9 bylFields et all (|2002l ) and ''U/'^U ~ 
1.2 by INakamura fc Shigevamal ( 2004 ) . by number. The ob- 
served ratio between the excess of ^Li and the ^Li abundance 
is '^Li/^Li ^ 15. We consider this ratio an upper limit since, 
given its fragility, some ^Li has probably been destroyed in 
previous evolutionary phases. The production of Be with- 
out a corresponding production of ®Li would be extremely 
difficult to understand. 

Nucleosynthetic calculations by INakamura et all (|200ll ) 
find the ejecta of HNe to have smaller amounts of C and O 
and large r amounts of Si, S , and Ar, w hen compared t o nor - 
mal SNe. INakamura et al.l ([200 ll ) and lNomoto et all (|2006l ) 
also note larger [(Zn,Co)/Fe] and smaller [(Mn,Cr)/Fe] ra- 
tios. An overabundance of Zn, which is not observed, can be 



avoided with a d eeper mas cuin, whi ch would also result in 
a larger [Ni/Fe] ([Nomoto et al.ll2006h . These are in qualita- 
tively agreement with the observations, supporting the HNe 
hypothesis. 

The weak s-process in massive stars seems to efficiently 
produce only elements w ith a mass numbeiQ up to 90 
([Ravet fc Hashimotd[2000l ). The Y overabundance may re- 
quire an enhanced flux of neutrons, which would also con- 
tribute to the production of Ni. The Ba overabundance, how- 
ever, is more difficult to understand. It is not clear whether 
this same mechanism would result in the overproduction of 
Ba. Moreover, a significant amount of Ba is expected to be 
produced only by the main s-process in AGBs or by the r- 
process, usually associated to massive stars. Since pollution 
by AGBs is not possible (see below), the Ba overabundance 
is likely a prod uct of a massive star. For example, although 
not expected, ;Mazzali,_Lucy fc Butler (1992) found Ba to be 
overabundant by a factor of 5 in the spectrum of SN 1987A. 
Although Ba might pose a problem for our scenario, we re- 
call that theoretical predictions for the r-process elements 
in HNe are not available. Therefore, whether this scenario is 
able to explain the Ba abundance is still an open question. 

The HNe scenario, at least qualitatively, is able to ex- 
plain most features observed in HD 106038 within a single 
peculiar event. More work, however, is still needed to show 
whether the scenario still holds quantitatively. A detailed 
comparison with nucleosynthetic predictions of theoretical 
models is necessary to validate or not the HNe hypothesis. 



3.2 Other scenarios 

In this subsection we present some alternative scenarios for 
the origin of the Be enhancement in HD 106038, which were 
discarded for the reasons presented below. 

A pollution by AGB stars, or any kind of evolved star, 
can at once be discarded. Although these may explain the 
s-process elements, and maybe Li, they do not explain the 
Be overabundance. On the contrary, the material ejected by 
an AGB or by a massive star, after successive mixing events, 
would be depleted in Be^ 

The SBs scenario ([Parizot|[2000l ) would be able to re- 
produce the observed Li/Be and Be/O ratios only with an 
extreme model where particles are accelerated from pure 
SNe ejecta. The SBs evolution models, however, predict that 
most material inside of a SB comes from the ISM. Moreover, 
the remaining chemical peculiarities are not typical of SNe 
II ejecta and thus can not be explained within the same 
scenario. 

Another possibility is the engulfing of a sub-stellar ob- 
ject, a planet or planetesim als debris as in star J37 of the 
open cluster NGC 6633 (Ash weU et"alll2005h . This, however, 
can be excluded with a robust quantitative argument. The 
accreted material would be confined to the surface convec- 
tive layer of the star. In a metal poor star this layer is much 
shallower tha n in a sola r metallicity star. With the equation 
given by Mur ray et al.l ([2001) we estimate the surface con- 
vective layer of HD 106038 to have ~ 4.5 x 10"^ Mq. The 



^ The coordinate, in mass, separating the part of the star that is 
ejected and the one that forms the remnant. 



^ The mass number of Y is A = 89 and of Ba is A 



136. 
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mass of Be in this layer is ~ 7.7 x 10~^^ M© while in a star 
with normal abundance of Be it would be ~ 4.8 x 10^^'' Mq. 
Assuming the accreted materijil to have a composition sim- 
ilar to chondrites meteorites ((Lodd crs 2003) a mass of Fe of 
5.3 X 10~^ Mq would also be accreted with the required mass 
of Be. However, the total mass of Fe in the convective layer 
of the star is ~ 3.3 x 10~^ M©. In this scenario all, or almost 
all, Fe in the convective layer of the star would come from 
the accreted material. HD 106038 would then originally be 
a population III metal free star; an extremely unlikely pos- 
sibility. In addition, we note that the most metal poor star 
found to host planets has [Fe/H] = —0.68 (jCochran et al.l 
I2OO7I ). 

A longer exposure to EPs would be the natural expla- 
nation if HD 106038, for some reason, was younger than halo 
stars of the same metallicity. Its Be abundance would be a 
result of the accumulated action of EPs in a cloud where 
star formation was, somehow, delayed. Its Be abundance 
higher than the solar photospheric one could be a sign of 
a solar or young e r age , in agreement with the suggestion 
of IPasguini et al. (|2004l ) that Be abundances could be used 
as a cosmic clock. The position of the star in the HR dia- 
gram, although not favourable for a good age determination, 
favours an older age and argues against this hypothesis. 

If the star originated in or near the Galactic centre, the 
enhanced star formation and supernovae events could pro- 
vide an enhanced EPs flux. This flux might also originate 
from a non-stellar source such as the central black hole. A 
bulge origin for the star, however, seems unlikely. In particu- 
lar the abundances of Ni, Y, and Ba are not compa tible with 
the ones of bulge stars (jMcWiUiam fc RichlflQOi ). Since it 
requires another source for the Ni and s-process overabun- 
dances, we also discard this hypothesis. 
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